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EXECUTIVE  SUMMARY 


The  Navy  Standard  Surf  Model,  SURF  3.0,  is  known  to  have  shortcomings  in  longshore  current  predic¬ 
tion.  Currents  are  often  found  to  be  too  weak  over  most  of  the  surf  zone  and  too  high  in  very  shallow  water. 
Both  bottom  friction  and  horizontal  eddy  diffusion  formulations  in  SURF  3.0  are  carefully  examined  and 
modified.  Instead  of  using  a  constant  friction  throughout  the  surf  zone,  a  new  bottom  friction  function  is 
developed  based  on  water  depth  and  surf  zone  width.  Data  acquired  during  the  1990  Duck  Experiment  on 
Low-frequency  and  Incident-band  Longshore  and  Across-shore  Hydrodynamics  (DELILAH)  and  the  1994 
DUCK94  experiment  are  used  to  derive  the  optimal  empirical  constants  for  the  bottom  friction  and  eddy 
diffusion  functions.  The  implementation  of  a  new  bottom-friction  function  and  eddy  viscosity  coefficient 
significantly  improves  the  magnitude  and  location  of  the  maximum  current,  and  greatly  reduces  the  unreal¬ 
istic  high  current  near  the  water  edge.  The  new  model  also  produces  lower  root-mean-square  errors.  The 
general  accuracy  of  resulting  longshore  current  estimates  is  further  verified  using  the  1980  data  set  from 
National  Sediment  Transport  Experiment  (NSTS),  conducted  at  Leadbetter  Beach,  Santa  Barbara,  Califor¬ 
nia.  The  improved  model  constitutes  significant  progress  in  longshore  current  prediction  and  will  be  used  to 
upgrade  the  SURF  3.0  model. 
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IMPROVEMENT  AND  VALIDATION  OF  THE 
NAVY  LONGSHORE  CURRENT  MODEL 


1.0  INTRODUCTION 

SURF  3.0  is  the  latest  version  of  the  Navy  Standard  Surf  Model  (NSSM).  It  is  a  streamlined  version  of 
SURF96  that  has  been  modified  by  Migues  et  al.  (1999)  to  conform  to  the  programming  standards  of  the 
Oceanographic  and  Meteorological  Master  Library  (OAML).  SURF3.0  was  delivered  to  the  System  Integra¬ 
tion  Department  of  the  U.S.  Naval  Oceanographic  Office  (NAVOCEANO)  in  February  1999.  SURF  3.0  has 
been  approved  by  NAVOCEANO  to  replace  SURF  2.1  as  the  operational  version  of  the  NSSM.  This  report 
presents  straightforward  yet  highly  effective  improvements  to  the  longshore  current  computations  in  SURF 
3.0. 


This  report  is  the  latest  in  a  recent  series  of  reports  that  address  various  aspects  of  the  NSSM,  including 
validation  and  alternative  initialization  methods. 

Mettlach  and  May  (1997)  examined  the  accuracy  of  model-derived  modified  surf  index  (MSI).  The 
model  used  (Mettlach  et  al.,  1996)  was  a  research  version  of  SURF  2.1.  It  was  shown  that  longshore  current 
is  a  relatively  high  contribution  to  the  MSI  and  that  errors  in  current  produce  concomitantly  high  errors  in 
MSI. 

May  and  Mettlach  (1997)  explored  the  utility  of  wave  model  estimates  and  depth  profile  estimates  using 
remote  sensing  techniques  for  SURF96  initialization.  The  accuracy  of  Navy  Wave  Model  (WAM)  nowcasts 
and  forecasts  for  surf  model  initialization  was  examined  in  detail. 

Hsu  et  al.  (1997)  conducted  a  thorough  evaluation  of  SURF96.  SURF96  and  SURF3.0  produce  nearly 
identical  results  when  certain  options  in  SURF96  conform  to  the  configuration  of  run  parameters  set  in 
SURF  3.0.  Their  model  was  initialized  and  evaluated  using  in  situ  meteorological,  oceanographic,  and  bathym¬ 
etry  data  acquired  during  two  experiments.  Four  cases  from  the  National  Sediment  Transport  Study  (NSTS) 
conducted  at  Leadbetter  Beach,  Santa  Barbara,  California,  during  February,  1980 — the  same  four  used  in 
this  study — were  used  to  examine  model  accuracy  with  a  planar  beach.  Seventy-nine  cases  from  the  Duck 
Experiment  on  Low-frequency  and  Incident-band  Longshore  and  Across-shore  Hydrodynamics  (PF.T.TT.AH), 
conducted  at  the  U.S.  Army  Corps  of  Engineers,  Field  Research  Facility  (FRF),  Duck,  North  Carolina, 
during  October  1990  were  used  to  examine  model  accuracy  with  a  barred  beach.  These  cases  are  also  used  in 
the  present  report.  Two  characteristic  shortcomings  in  longshore  current  estimates  were  discovered: 

•  Over  most  of  the  surf  zone,  the  calculated  longshore  current  was  approximately  half  of  that  mea¬ 
sured. 

•  Very  near  the  edge  of  the  water,  in  water  less  than  1  -m  depth,  the  model  often  produced  a  second  peak 
that  was  larger  than  the  first  peak,  contrary  to  measured  data.  This  problem  is  particularly  serious 
over  a  planar  beach.  In  the  present  report,  methods  to  correct  both  of  these  shortcomings  are  given. 
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Based  on  SURF96,  Mettlach  et  al.  (1999)  estimated  surf  zone  width.  More  than  600  surf  zone  width 
determinations  from  10-min  averaged  video  records  were  optimally  correlated  to  the  most  offshore  location 
where  the  surf  model  estimated  10  percent  wave  breaking  to  occur.  This  simple  method  produced  mean  surf 
zone  width  errors  of -2  m,  with  a  standard  deviation  of  30  m.  Video-determined  widths  ranged  from  5  to  408 
m.  In  the  present  study,  it  is  shown  that  model-derived  surf  zone  width  can  be  effectively  used  to  improve 
longshore  current  estimates  in  SURF3.0. 

Lundberg  et  al.  (1997,  1999)  have  developed  a  semi-empirical  longshore  current  model  that  does  not 
need  bottom  depth  as  input.  In  their  model,  the  shape  of  the  longshore  current  profile  in  the  cross-shore 
direction  is  assumed.  The  current  is  near  zero  at  a  distance  of  two  surf  zone  widths  from  the  beach,  steadily 
increases  to  a  maximum  at  a  distance  of  one-half  surf  zone  width,  and  then  steeply  decreases  to  zero  at  the 
beach.  Model  skill  is  derived  largely  from  the  fact  that  longshore  current  profiles  generally  match  this  bell 
shape  over  both  planar  beaches  and  barred  beaches. 

This  report  shows  that  longshore  current  estimates  from  SURF  3.0  can  be  improved  considerably  by 
replacing  a  constant  bottom  friction  coefficient  with  a  variable  bottom  friction  function  based  on  surf  zone 
width.  The  empirical  constants  for  bottom  friction  and  horizontal  eddy  viscosity  are  derived  based  on  field 
data.  The  improved  model  is  proven  to  be  robust. 

2.0  DESCRIPTION  OF  SURF  3.0 

This  section  describes  wave  transformation,  wave  setup,  and  longshore  current  formulations  used  in 
SURF  3.0.  The  main  assumption  governing  these  equations  is  that  there  is  along-shore  uniformity  in  depth. 
Numerical  methods  for  computing  the  solutions  of  the  differential  equations  representing  wave  transforma¬ 
tion  and  longshore  current  are  described  in  detail.  Empirical  constants  used  in  the  model  are  given. 

Wave  breaking  and  transformation  is  based  on  the  concept  that  organized  wave  motion  at  breaking  is 
converted  to  a  turbulent  roller  riding  over  the  underlying  wave.  The  production  and  advection  of  wave 
energy  balances  the  dissipation  of  wave  energy.  The  equations  for  wave  transformation  using  roller  theory 
are  taken  from  Lippmann  et  al.  (1996),  which  is  based  on  the  parametric  wave  transformation  model  of 
Thornton  and  Guza  (1983). 

The  longshore  current  model  is  based  on  the  radiation  stress  theory  of  Longuet-Higgins  (1970a,  1 970b), 
with  an  added  term  for  wind  stress  given  in  Earle  (1991).  The  numerical  solution  for  longshore  current  uses 
Gaussian  elimination  of  a  tri-diagonal  matrix,  which  is  given  in  Kraus  and  Larson  (1991). 

2.1  Wave  Transformation 

The  energy  flux  equation  balances  the  wave  and  roller  fluxes  of  energy  density  to  energy  dissipation  £, 

d(EwCg  cosO)  ^  d(Erc cosQ)  _  ^  (1) 

dx  dx 

The  energy  of  incoming  waves  travels  at  the  wave  group  velocity  cg,  and  rollers,  which  are  associated 
with  depth-limited  breaking,  are  assumed  to  travel  at  the  phase  velocity  of  the  wave  c.  Using  linear  wave 
theory,  the  wave  energy  per  unit  area  of  sea  surface  in  a  distribution  of  waves  with  root-mean-square  wave 
height  H  is  given  as 
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g  P S^rms  ’  (2) 

with  g  being  the  acceleration  due  to  gravity  (9.81  m/s2)  and  p  the  density  of  seawater  (1030  kg/m3).  Linear 
wave  theory  also  yields  the  following  relationships  among  wave  energy  speed  c  ,  angular  frequency  co, 
radian  wave  number  k,  wave  phase  speed  c,  depth  h,  wave  frequency  f  and  wavelength  L : 


_d(£>_  fl  kh  \ 

8  dk  V2  sinh  2khf 

(3) 

CO 

c  —  —  =  fL; 
k 

(4) 

oo2  =  gk  tanh  kh . 

(5) 

Using  the  assumption  that  depth  contours  in  the  nearshore  are  locally  straight  and  parallel  to  the  beach, 
then  Snell’s  law  can  be  used  to  compute  wave  direction  0  if  the  offshore  wave  angle  0o  and  its  frequency  and 
depth  are  known. 


sin#  _  sin #0 

c  c0 


(6) 


The  surf  model  is  initialized  with  the  measured  or  estimated  wave  height,  period,  and  direction  at  an 
offshore  point.  These  wave  parameters  are  best  derived  from  a  complete  directional  wave  spectrum  S(f,Q). 
The  root-mean-square  wave  height  H rms  is  computed  from  the  zero  moment  of  wave  variance. 


Hrms  ~ 


LJ 

11  mo  _ 

V2  " 


71/2 
8  1 
-71/2 


^2 


J  S(f,Q)dfdQ 
o 


The  peak  wave  frequency  in  the  spectrum  is  that  frequency  band  with  a  maximum  in  energy. 


(7) 


f 

’  7t/2  T 

fP=f 

max 

J  S(/,0>/0 

V 

_  7l/2 

(8) 


It  is  used  for  model  initialization  and  is  assumed  conserved  throughout  the  surf  zone. 

The  wave  angle  used  for  model  initialization  0o  is  computed  from  the  total  energy  of  the  input  wave 
spectrum  and  the  total  energy  of  the  along  shore  component  of  cross  shore  radiation  stress  S  i.e., 


A  1  •  -1 

#0  =-sm 


k/2 


2  j  js^mdfde 


-nil  0 


nl 2 

j  \S(f,0)dfd9 


-nil  0 


1  + 


2k(fp,h)h 
sinh2  k(fp,h)h 


1\ 


J/ 
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Jxy 


(/,0)  =  S(/,0) 


1 

-  + 
2 


k(f,h)h 

sinh2  k(f,h)h 


cos  0  sin  0. 


(10) 


Because  of  wave  refraction,  most  wave  angles  will  be  within  45  deg  of  shore  normal.  Because  of  the 
relatively  small  wave  angles  in  the  surf  zone,  errors  in  wave  angle  produce  relatively  greater  errors  in  longshore 
current  than  in  wave  height  because  wave  height  is  proportional  to  the  cosine  of  the  wave  angle  while  the 
longshore  current  is  proportional  to  its  sine.  If,  for  example,  the  true  wave  angle  is  20  deg  but  the  wave  angle 
used  in  the  model  is  10  deg,  then  there  will  be  a  4.7  percent  error  in  wave  height,  but  a  -49.2  percent  error  in 
longshore  current.  It  is  stressed  that  changes  in  the  definition  of  wave  angle  may  have  significant  effects  on 
longshore  current  estimates. 

Surface  roller  energy  density  is  given  as 


Er=-prc2-.  (11) 

r  2 Kr  L 

This  equation  is  derived  from  the  kinetic  energy  density  of  the  roller  volume,  i.e.,  from  mv2/ 2,  where  v  is 
assumed  to  be  the  speed  of  travel  of  the  roller,  which  is  the  phase  speed  of  the  underlying  wave.  The  kinetic 
energy  density  is  found  relative  to  a  unit  wavelength  and  the  surface  area  A  of  the  roller  covering  the  face  of 
the  underlying,  breaking  wave, 


4/ztanct 


(12) 


This  area  is  approximated  from  the  root-mean-square  wave  height  at  the  time  of  wave  breaking  Hb,  local 
bottom  depth  h,  and  the  angle  of  the  wave-roller  stress  vector  o,  which  gives  the  slope  of  the  wave  front.  The 
density  of  aerated  water  comprising  the  roller  pr  is  less  than  that  of  the  water  in  the  underlying  wave  below 
the  roller,  although  the  density  of  seawater  p  is  used  as  a  simplification.  Hb  is  actually  Hrms  in  model  calcula¬ 
tions  of  A. 

Energy  dissipation  is  assumed  to  occur  exclusively  at  the  shear  surface  interface  between  the  roller  and 
the  organized,  underlying  wave.  Bottom  friction  is  neglected.  The  turbulent  kinetic  energy  (TKE)  per  unit 
area  produced  at  the  wave-roller  interface  is  given  as 

xi=pgAsino.  (13) 

Thus,  energy  loss  per  horizontal  surface  area  and  per  unit  time  e  (with  units  W/m2)  is  given  as 


3 

£  =  ^-  =  —  pgf„—  COSC.  (14) 

L  4  p  h 

Energy  dissipation  is  calculated  by  accounting  for  changes  in  breaker  wave  height  Hb  (=  Hrm)  throughout  the 
surf  zone,  which  are  functionally  related  to  the  probability  density  that  a  wave  of  a  given  height  is  either 
breaking  or  broken.  The  probability  density  pb(H )  that  a  wave  of  given  height  H  is  either  breaking  or  broken. 
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pb(H)  =  W(H)p(H), 

is  found  from  a  Rayleigh  probability  density  distribution  p(H),  i.e., 


(15) 


H 

p(H)  =  -^-e  1 HrmsJ 

TJ / 

Ll  f.yy-  r 


(16) 


multiplied  by  a  factor  W(H),  which  is  an  empirically  derived  weighting  function.  In  SURF  3.0  the  following 
weighting  function  is  used, 


( H  \4 

W(H)= 

l  yh  J 


(  (H^ 


1-e 


yhj 


<i, 


(17) 


where  y  =  0.42  represents  a  constant  ratio  between  and  h.  This  weighting  function  is  taken  from  Thornton 

and  Guza  (1983),  who  derived  it  from  a  limited  set  of  wave  distributions  that  contained  breaking  and 
nonbreaking  waves. 

In  the  model,  the  cubed  breaker  wave  height  is  calculated  numerically  using  a  trapezoidal  integration 
from  0  to  5 H  , 

rms’ 

oo  5// 

H3b  =  lH3pb(H)dH=  X  H3pb(H)AH.  (18) 

0  H= 0 

The  energy  flux  balance  equation  is  solved  using  a  forward-stepping  numerical  function  from  offshore  to  the 
edge  of  the  water  in  increments  of  Ax , 

(Ewcg  cos  6)2  +(£rccos0)2  =  Axe  +  (£wcg  cos0)j  +(£rccos0)! .  (19) 

Subscripts  1  and  2  indicate  points  in  cross-shore  direction,  with  1  being  farther  offshore  than  2.  The  equation 
is  balanced  by  finding  Hrms2  with  an  iterative  solution. 

2.2  Wave  Setup 

Setup  of  the  mean  water  level  r\  is  given  by 


d$xx,w  j 
dx 


d$xx,r 

dx 


+Pg/3=o. 

dx 


(20) 


The  first  term  is  the  gradient  in  radiation  stress  due  to  wave  motion,  and  the  second  term  is  that  due  to  roller 
motion.  The  third  term  is  the  pressure  gradient  from  setup.  The  wave  contribution  to  radiation  stress  is 


(21) 


6 


Hsu,  Mettlach,  and  Earle 


The  roller  momentum  flux  is 


=2£rcos20.  (22) 

In  SURF  3.0,  the  setup  is  added  to  the  mean  water  level  to  compute  wave  height.  A  simple  forward  finite- 
difference  numerical  function  is  used  to  compute  setup. 

2.3  Longshore  Current 

The  longshore  current  model  used  in  SURF  3.0  balances  the  longshore  (y-directed)  momentum  flux 
from  waves  in  the  surf  zone  %H,  the  horizontal  exchange  of  momentum  due  to  turbulent  mixing,  the  dissipa¬ 
tion  of  momentum  due  to  bottom  friction  xfl,  and  the  generation  of  momentum  from  the  wind  x(;. 


d 

T"+p* 


dV\ 

h\t—  \-XB  +  Ty  -0. 
dxj 


(23) 


Wave  input  for  the  momentum  balance  equation  is  given  by 

sinO 
xH  =e - . 


(24) 


In  the  second  term  of  the  momentum  balance  equation,  lateral  or  horizontal  eddy  viscosity  is  represented  by 
(Battjes,  1975) 


v,  =  Mh\ 


^V/3 

U>, 


such  that  Mis  an  empirical  constant,  and  is  set  to  2  in  SURF  3.0. 
The  bottom  shear  stress  term  is  modeled  simply  as 


(25) 


xB-CfpuV.  (26) 

The  bottom  drag  coefficient  cf  =  0.007  is  constant  throughout  the  surf  zone.  It  is  shown  later  that  using  a 
variable  bottom  friction  coefficient  function  improves  model  estimates  compared  to  field  data.  This  is  intu¬ 
itively  reasonable  in  that  there  is  often  preferential  sorting  of  sediment  material  by  grain  size  in  the  surf 
zone,  yielding  areas  of  relatively  coarse  sand  in  the  swash  zone  and  progressively  finer  sand  with  greater 
depth  and  distance  offshore. 

The  variable  u  is  the  magnitude  of  the  near-bottom  horizontal  wave  orbital  velocity, 


rms 


U  = 


sinh  kh 


(27) 
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Finally,  the  wind  stress  from  a  wind  of  speed  U  and  direction  0  is 


T u  —  U  sin  Ovpacd, 


pa=  1.2  g/m3,  and 


1.2875  xlO~3  U<  7.5m /s 

cd=  ,  .  (29) 

(0.8  +  0.06517)1 0~3  U>  7.5  m /s 

The  longshore  current  is  computed  numerically  using  three  passes  across  the  cross-shore  grid  represent¬ 
ing  the  surf  zone.  In  the  first  pass,  energy  dissipation  and  wave  height  are  computed.  Three  intermediate 
variables  A.,  B  ,  and  Ci  are  also  computed  during  this  pass,  which  proceeds  from  offshore,  where  i  =1,  to  near 
the  water’s  edge, "There  i  -  N.  These  intermediate  variables  are  given  as: 


A,-  =  Mphf 


cfPfp^ftHrms,i 


sinh  kjhj 


_  -E.  sine,  _u2 


-U  sin  0[/pacrf. 


In  the  second  pass,  the  two  intermediate  variables  E.  and  F.  are  computed.  This  pass  proceeds  from  near 
the  beach,  where  i  =  N  -1  to  offshore,  where  i  =  2: 


b,ax2  +  a.,+A(i-£„,)’ 


p  _ _ AiFi+\  CjAx  _ 

'  BjAx2  +  Aj_i  +  Ai  (l  -  Ei+i ) ' 


At  the  beach, 


Em  —  0 


r  _  £n  sin60 

b N  - 

BN  co 


(36) 
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The  third  pass  yields  the  model’s  estimate  of  longshore  current  V.  from  i  =  2  to  /  =  N.  This  pass  proceeds 
from  offshore  to  near  the  beach, 


Vi=ElVl_l+Fi. 

This  pass  begins  with  the  steady-state  wind-generated  current  Vv 

Vj  =  f—  —UsinOy . 
VP  cf 


(37) 


(38) 


3.0  MODEL  EVALUATION  OF  SURF  3.0 

This  section  summarizes  an  evaluation  of  the  accuracy  of  SURF  3.0  in  estimating  H  and  V.  This  is  done 
using  wave  height  and  longshore  current  measurements  acquired  during  DELILAH  and  DUCK94.  The 
model  is  initialized  with  depth,  tide,  and  offshore  wave  and  wind  measurements  from  these  experiments. 

Because  of  the  demands  that  are  placed  on  the  operational  performance  of  the  NSSM,  it  is  imperative 
that  it  be  tested  under  a  wide  range  of  wind  and  wave  conditions.  The  combined  DELILAH-DUCK94  data 
set  is  used  for  such  testing.  Of  the  cases  selected  from  these  two  field  experiments,  the  wind  speed  ranged 
from  1  to  12  m/s,  and  the  offshore  significant  wave  height  ranged  from  0.4  to  4.0  m.  Figures  1  and  2  give 
time  series  of  wind  speed,  wave  height,  wave  direction,  and  peak  wave  period  of  the  cases  used  in  this  study. 


Pig.  i  — Wind  and  wave  measurements  for  79  DELILAH  cases — wind  speed,  significant 
offshore  wave  height,  peak  wave  direction,  and  peak  wave  period 
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Fig.  2  —  Wind  and  wave  measurements  for  144  DUCK94  cases — wind  speed,  significant 
offshore  wave  height,  peak  wave  direction,  and  peak  wave  period 


The  beach  at  the  FRF  is  representative  of  many  beaches  of  interest  to  the  Navy.  The  nearshore  gradient 
number  at  Duck,  depending  on  the  along-shore  location,  ranges  from  approximately  85  to  105,  with  a  mean 
of  94.  The  gradient  number  is  a  single  dimensionless  number  Ax/30,  where  Ax  is  the  distance  from  the  beach 
in  feet  to  a  depth  of  30  ft.  Using  an  unclassified  listing  based  on  the  Defense  Planning  Guide  (DPG),  which 
contains  distance  to  depth  of  30  ft  for  248  potential  landing  sites,  Mettlach  et  al.  (1999)  have  shown  that  the 
median  beach  gradient  number  in  the  DPG  data  set  is  92,  in  which  the  values  range  from  8  to  241 1 .  Approxi¬ 
mately  70%  of  the  DPG  gradient  numbers  are  between  25  and  200.  This  median  value  is  very  close  to  that  at 
Duck. 

The  conditions  during  the  two  experiments  are  markedly  different  in  one  particular  way — the  averaged 
peak  wave  direction  0p  during  DELILAH  was  from  —24.9  deg  (positive  counterclockwise  from  shore  nor¬ 
mal,  seaward  (070  deg  north),  while  during  DUCK94,  the  averaged  peak  wave  direction  was  +2.6  deg.  The 
peak  wave  direction  is  defined  as  that  direction  in  which  there  is  a  maximum  in  wave  energy.  During 
DELILAH,  the  long  pier  south  of  the  surf  zone  observation  area  may  have  affected  longshore  current  devel¬ 
opment  from  these  waves  arriving  from  the  southeast.  During  DUCK94,  the  frequent  occurrence  of  waves 
directly  onshore,  of  course,  produces  very  small  longshore  currents.  Also,  very  high  winds  during  DUCK94 
could  have  set  up  non-wave-driven  strong  coastal  currents  in  ways  that  are  not  considered  by  the  model 
physics. 

3.1  Wave  Height 

Figure  3  is  a  scatter  plot  of  the  root-mean-square  wave  heights  acquired  during  DELILAH  and  DUCK94 
vs  the  corresponding  model  estimates.  The  sensor  description  and  arrangement  is  easily  accessible  from  the 
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Pig,  3  —  Scatter  diagram  of  measured  root-mean-square  wave  heights  vs  estimates  of  wave 
height.  Dots  denote  DELILAH  measurements;  circles  denote  DUCK94  measurements.  N  = 
1796  observations.  Skill  statistics  based  on  model  minus  measured  differences  of  set  from 
combined  data  sets  are:  linear  correlation  coefficient  R  =  0.945;  minimum  error  emfa  =  -1.0 1  m; 
mean  error  E  =  0.04  m;  maximum  error  e  =  0.71  m;  standard  deviation  of  the  errors  a  = 

mean  wax 

0.15  m;  mean  fractional  error  e100  =  6.2%;  mean  of  the  absolute  values  of  the  fractional  error 
e  ^ .  =  12.5  %.  The  solid  line  is  a  linear  regression  curve  with  slope  m  =  0.90,  and  ^-intercept 

lOOabs 

b  =  0.15  m. 


Corps  of  Engineers  Web  site,  therefore  is  not  repeated  here.  The  skill  of  the  model  is  quite  high,  with  a  mean 
percent  error  of  +6.2  percent,  which  is  close  to  the  minimum  error  found  by  Lippmann  et  al.  (1996)  using 
just  nine  cases.  The  model  was  run  with  the  standard  model  parameters  given  in  Section  3.0,  i.e.,  y  =  0.42 
and  a  =  5°.  For  this  validation,  the  Self  Start  option  of  the  model  (see  Hsu  et  al.,  1997)  was  used.  The  Self 
Start  option  suppresses  wave  energy  dissipation  until  a  threshold  of  5  percent  wave  breaking  occurs.  It  is 
noted  that  the  slope  of  the  best-fitted  linear  regression  curve  through  the  measurement-model  pairs  has  a  y- 
intercept  of +0.15  m/s  and  a  slope  m  of  0.90,  which  suggests  that  true  wave  heights  above  1.5  m  are  under¬ 
estimated  by  the  model. 

An  additional  series  of  15  sets  of  surf  model  runs  using  the  DELILAH  and  DUCK94  data  sets  were 
made  to  examine  the  sensitivity  of  the  model  to  its  two  wave  energy  dissipation  parameters — y  and  o  [Eq. 
(17)  and  Eqs.  (12)-(14)],  respectively).  This  examination  is  analogous  to  that  made  by  Lippmann  et  al. 
(1996),  who  calculated  the  percentage  of  error  between  modeled  and  measured  wave  height  for  individual 
surf  zone  measurement  cases.  It  the  present  examination  the  percentage  of  error  and  several  other  skill 
statistics  are  calculated  for  combinations  of  y  and  o  for  y  ranging  from  0.3  to  0.5  and  a  ranging  from  2°  to  6°. 
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The  results  are  given  in  Table  1 ,  which  shows  that  the  optimal  combination  of  y  and  a  is  y=  0.40  and  o  =  6°, 
which  is  very  close  to  the  standard  values  used  in  the  model.  Note  that  unlike  the  set  of  runs  represented  in 
Fig.  3,  the  additional  sensitivity  study  has  been  done  without  the  Self  Start  option,  allowing  a  uniform 
number  of  validating  measurements  for  each  set  of  runs  ( N  =  2660).  The  higher  skill  for  y  =  0.42  and  c  =  5° 
in  Table  1  (than  given  in  the  caption  of  Fig.  3)  is  due  to  the  greater  number  of  wave  measurements  farther 
offshore  being  incorporated  into  the  validation  (than  the  fewer  number  of  measurements  (N =1796)  from  the 
runs  with  the  Self  Start  option). 

3.2  Longshore  Current 


Figure  4  is  a  similar  plot  of  longshore  current  measurements  from  DELILAH  and  DUCK94  vs  the 
corresponding  model  estimate.  The  skill  of  the  model  in  estimating  V  is  much  less  than  that  of  the  model  in 
estimating  H rms-  The  most  striking  feature  is  the  low  value  of  slope  (m  =  0.48),  indicating  underprediction  of 
currents.  The  statistics  listed  under  the  figure  caption  are  discussed  in  Section  5.2. 


Table  1  —  Skill  Statistics  from  15  Sets  of  Surf  Model  Run 
Using  Values  of  y  and  o  (°)  Noted. 


Y 

a 

N 

R 

£ 

£ 

£ 

a  £  £ 

m 

b 

mm 

mean 

max 

100 

j.00abs 

30 

2 

2660 

0.942 

-1.038 

-0.112 

0.614 

0.165  -6.832 

13.36 

0.768 

0.131 

30 

4 

2660 

0.946 

-1.038 

-0.130 

0.584 

0.170  -8.212 

14.35 

0.723 

0.160 

30 

6 

2660 

0.946 

-1.038 

-0.136 

0.553 

0.175  -8.703 

14.87 

0.707 

0.170 

36 

5 

2660 

0.960 

-1.007 

-0.045 

0.645 

0.137  -0.663 

10.46 

0.821 

0.143 

37 

5 

2660 

0.961 

-1.007 

-0.042 

0.645 

0.137  -0.381 

10.44 

0.822 

0.145 

38 

5 

2660 

0.961 

-1.007 

-0.039 

0.675 

0.136  -0.094 

10.48 

0.822 

0.147 

39 

5 

2660 

0.961 

-1.007 

-0.037 

0.675 

0.136  0.186 

10.53 

0.823 

0.148 

40 

2 

2660 

0.951 

-1.007 

0.018 

0.706 

0.141  4.692 

11.56 

0.908 

0.113 

40 

5 

2660 

0.961 

-1.007 

-0.034 

0.675 

0.136  0.437 

10.59 

0.824 

0.150 

40 

6 

2660 

0.962 

-1.007 

0.006 

0.675 

0.126  3.753 

10.27 

0.882 

0.130 

41 

5 

2660 

0.961 

-1.007 

-0.032 

0.706 

0.136  0.714 

10.70 

0.824 

0.152 

42 

5 

2660 

0.960 

-1.007 

-0.030 

0.706 

0.136  0.958 

10.81 

0.825 

0.154 

50 

2 

2660 

0.942 

-1.007 

0.117 

0.797 

0.166  13.523 

16.56 

1.017 

0.099 

50 

4 

2660 

0.951 

-1.007 

0.115 

0.767 

0.152  13.257 

16.08 

1.014 

0.100 

50 

6 

2660 

0.953 

-1.007 

0.114 

0.767 

0.148  13.120 

15.86 

1.014 

0.099 

y  =  coefficient  used  with  Thornton  and  Guza  wave  dissipation  model 
o  =  angle  of  the  wave-roller  stress  vector  (degrees) 

N  =  number  if  model-measurements  pairs 

R  =  linear  correlation  coefficient  between  model  estimates  and  measurements 

e  =  minimum  error  (m) 

min  v  y 

£ mean  =  simple  mean  of  the  individual  errors  (model  -  measurement)  (m) 

e  =  maximum  error  (m) 

a  =  standard  deviation  of  the  error  (±  m) 

£100  =  100  x  mean  fractional  error  ([model -meas.])/meas.;  meas.  *  0) 

eiooabs  =  100  x  mean  absolute  value  of  fractional  error  ([model-meas.]/meas.;  meas.  *=  0) 

m  =  slope  of  regression  line  through  model-measurement  pairs 

b  =  ^-intercept.  Offset  of  regression  line  through  model-measurement  pairs  (m) 


N  is  number  of  measurement-model  estimate  pairs  from  the  set  of  surf  model  runs.  Linear  correlation 
coefficient  is  R.  Minimum  error  is  Emjn  (m).  Mean  error  is  zmean  (m).  Maximum  error  is  ^(m).  Stan¬ 
dard  deviation  of  error  is  a  (m),  which  should  not  be  confused  with  the  angle  c.  Mean  fractional  error 
in  percent  is  £100.  Mean  of  the  absolute  values  of  the  fractional  error  in  percent  is  £100abs  (%)  Slope  of 
regression  curve  through  model-measurement  pairs  is  m.  Y-intercept  of  regression  curve  is  b  (m). 
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c  -  0 . 0070;  M=2 . 00 


Fig.  4  —  Scatter  diagram  of  measured  longshore  currents  vs  corresponding  model 
estimates  of  current.  Dots  denote  DELILAH  measurements;  circles  denote  DUCK94 
measurements.  N  =  1930  observations.  Skill  statistics  based  on  model  minus 
measured  differences  of  set  from  combined  data  sets  are:  linear  correlation  coefficient 
R  =  0.87;  minimum  error  e  =  -1.15  m/s;  mean  error  e  _  =  -0.10  m/s;  maximum 
error  e  =  0.57  m/s;  standard  deviation  of  the  errors  a  =  0.27  m/s;  mean  fractional 

max 

error  eioo  =  -4-9.9  %;  mean  of  the  absolute  values  of  the  fractional  error  E100abs=  97 .5 
%.  The  solid  line  is  a  linear  regression  curve  with  slope /n  =  0.48  and  y-intercept 
b  =  -0.07  m/s. 


Accurate  estimates  of  the  maximum  current  in  the  surf  zone  Vmax  are  important  with  regard  to  modified 
surf  index  (MSI)  computations.  MSI,  defined  in  the  Joint  Surf  Manual  (Commander,  Naval  Surface  Force, 
Pacific  and  Commander,  Naval  Surface  Force,  Atlantic,  1987)  is  the  sum  of  weighted,  nondimensional 
values  of  surf  parameters.  These  include  wave  height,  wave  period,  wave  direction  relative  to  the  beach,  Vmax, 
breaker  type  and  wind  speed.  Each  foot  of  wave  height  contributes  1  MSI  unit  to  the  total.  Each  knot  of 
longshore  current  contributes  3  MSI  units.  In  Fig.  5,  measured  V ^  from  the  DELILAH  and  DUCK94  cases 
are  plotted  against  the  corresponding  SURF  3.0-modeled  Vmax.  The  statistics  of  Fig.  5  are  similar  to  Fig.  4.  It 
should  be  noted  that  the  slightly  improved  values  are  due  to  the  fact  that  some  computed  V^are  located  at 
the  end  of  the  beach,  far  from  the  measured  peak  near  the  middle  of  surf  zone.  Those  values  probably  should 
not  be  considered  as  valid  values  in  computing  the  statistics. 

4.0  PARAMETERIZATION S  IN  LONGSHORE  CURRENT  MODEL 

Physical  processes  in  the  surf  zone  are  especially  complex  and  require  some  degree  of  parameterization 
for  efficient  modeling.  The  set  of  equations  given  in  Section  2.0  relies  on  several  empirical  parameters  and 
parameterizations  for  closure,  two  of  which  are  discussed  in  this  section:  c  and  v(. 
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Fig.  5  —  Scatter  diagram  of  maximum  measured  longshore  currents  vs  the 
corresponding  maximum  model  estimates  of  current.  Dots  denote  DELILAH 
measurements;  circles  denote  DUCK94  measurements.  N  =  223  observations.  Skill 
statistics  based  on  model  minus  measured  differences  of  set  from  combined  data  sets 
are:  linear  correlation  coefficient/?  =  0.86;  minimum  error  Bmjn  =  -0.93  m/s;  mean  error 
e  =  -0, 19  m/s;  maximum  error  e  =  0.48  m/s;  standard  deviation  of  the  errors  a  = 
0.32  m/s;  mean  fractional  error  ei00  =  -^40. 1  %;  mean  of  the  absolute  values  of  the 
fractional  error  e100abs  =  74.2%.  The  solid  line  is  a  linear  regression  curve  with  slope 
m  =  0.63  and  ^-intercept  b  -  -0.1  m/s 


4.1  Bottom  Friction  Coefficient  cf 

The  value  of  the  bottom  friction  coefficient  cf  is  the  main  fitting  parameter  in  longshore  current  model¬ 
ing.  The  value  used  in  SURF  3.0  is  the  same  value  used  in  previous  versions  of  the  navy  model,  i.e.,  cf  = 
0.007,  constant  throughout  the  surf  zone  (Earle,  1991). 

Thornton  and  Guza  (1986)  in  citing  Shemdin  et  al.  (1978)  state  that  reported  values  of  cf  range  between 
0.005  to  0.3.  Thornton  and  Guza,  under  the  assumption  of  a  constant  cyin  the  surf  zone,  determined  optimal 
values  of  c.from  four  field  cases  by  solving  through  trial  and  error  the  least-square  error  between  calculated 
and  measured  longshore  current  values  shoreward  of  the  mean  breaker  line.  For  a  linear  bottom  stress  for¬ 
mulation  of  longshore  current,  similar  to  that  used  in  SURF96,  their  optimal  values  ranged  from  0.005  to 
0.009.  For  the  same  four  cases,  using  a  similar  linear  bottom  stress  model,  Kraus  and  Larson  (1991)  found 
optimal  values  ranging  from  0.004  to  0.006.  Reniers  and  Battjes  (1997)  in  a  laboratory  study  of  longshore 
currents  over  a  barred  beach  found  that  cf=  0.007,  constant  was  an  optimal  value  for  a  linear  model  incorpo¬ 
rating  rollers  in  the  energy  dissipation  calculations. 
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Smith  et  al.  (1993)  examined  three  candidate  forms  of  cf  The  first  was  to  assume  that  cf  is  a  constant 
value.  The  second  was  to  use  the  wave  friction  coefficient  of  Jonsson  (1966),  which  was  used  by  Thornton 
(1970)  and  is  given  by 


1 


t 

—  — 0.08  4-  log — , 


(39) 


in  which  £  is  the  maximum  water  particle  excursion  amplitude  of  fluid  motion  at  the  bottom  as  predicted  by 
linear  wave  theory,  i.e., 


H  1 
2  sinh/r/z’ 


(40) 


and  r  is  a  roughness  parameter,  being  the  measure  of  ripple  height.  Third,  the  Manning  friction  coefficient 


c/  —  g—y  (41) 

dA 

was  used.  In  this  formulation,  n  is  an  empirically  determined  coefficient,  called  the  Manning  resistance 
coefficient,  and  d  is  the  hydraulic  radius  of  wide  open-channel  flow.  Of  the  three  formulations,  the  Manning 
friction  coefficient  was  found  to  produce  the  least  error  in  longshore  current.  Smith  et  al.  (1993)  do  not 
clearly  define  d ,  but  do  show  in  a  plot  from  DELILAH  (14  October  1990  0100  EST)  the  Manning  friction 
coefficient  ranging  from  cf  =  0.0022  far  offshore,  where  h  >  10  m  to  cf>  0.06  near  the  edge  of  the  water.  The 
lower  value  is  less  than  those  mentioned  in  Thornton  and  Guza  (1986)  or  Kraus  and  Larson  (1991). 

Sancho  (1999),  who  examined  unsteady  nearshore  currents  on  longshore  varying  topographies  with 
DELILAH  measurements,  used  an  approximation  of  Eq.  (39)  based  on  Swart  (1974).  For  the  0400  UTC  10 
October  1990  case,  estimates  for  pranged  from  0.006  offshore  to  0.017  at  the  shoreline. 

During  the  SUPERDUCK  experiment  conducted  at  the  FRF  in  October  1986,  Whitford  and  Thornton 
(1986)  derived  quantitative  estimates  of  various  terms  in  the  momentum  equation  governing  the  mean  along¬ 
shore  flow  in  the  surf  zone.  An  instrumented  sled  acquired  mean  current  and  bottom  shear  stress  measure¬ 
ments  at  various  locations  across  the  surf  zone.  Coefficients  of  c  were  estimated  to  range  from  0.004  sea¬ 
ward  of  the  bar,  to  0.003  at  the  crest  of  the  bar,  to  0.001  over  the  trough  shoreward  of  the  bar.  It  is  shown  that 
the  shoreward  decrease  in  cf  is  statistically  significant  and  suggested  to  have  been  due  to  decreasing  wave 
action  shoreward. 

More  recently,  Garcez  Faria  et  al.  (1998)  have  produced  another  formulation  for  c^from  field  data,  one 
that  is  dependent  on  depth  and  k  ,  an  apparent  bed  roughness  length  scale: 


c 


f 


(42) 


This  formula  is  based  on  field  measurements  of  ^acquired  at  Duck  during  three  strong  current  days  during 
DUCK94.  cf  ranged  by  more  than  an  order  of  magnitude  (0.0006  -  0.0120).  As  shown  in  Fig.  6,  cf  tends  to 
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Fig.  6  —  Measurements  of  ^acquired  by  Garcez  Faria  et  al.  (1998)  at  Duck,  North  Carolina, 
during  DUCK94.  (Cross-shore  distance  is  relative  to  a  fixed  location  at  the  Field  Research 
Facility  (FRF)  at  Duck,  North  Carolina.)  Depending  on  tide,  the  edge  of  the  water  is  located  at 
100-120  m. 

increase  shoreward  due  to  increasing  bottom  form  irregularities  shoreward.  There  was  a  noticeable  transi¬ 
tion  from  a  smooth,  near-planar  bottom  form  seaward  of  the  bar,  to  rough  mega  ripples  over  the  bar,  to  less 
rough  wave  ripples  very  near  the  beach.  Unfortunately,  Eq.  (42)  cannot  be  easily  used  in  numerical  models 
because  apparent  bed  roughness  cannot  be  known  without  detailed  velocity  profile  measurement. 

In  view  of  the  wide  range  of  formulations  for  bottom  friction  and  its  great  importance  in  longshore 
current  estimates,  which  is  shown  later  in  this  report,  a  practical  formulation  for  c  is  developed.  This  formu¬ 
lation  greatly  reduces  spuriously  high  estimates  of  V  in  very  shallow  water  and  increases  V  to  values  much 
closer  to  measurements  over  the  great  majority  of  the  surf  zone. 

4.2  Parameterization  of  Lateral  Eddy  Viscosity  Vr 

The  lateral  eddy  viscosity  V(  is  also  an  important  empirical  parameter  in  longshore  current  models. 
Thornton  and  Guza  (1986)  used 


v,  =  AM  JC I  Jgh  (43) 

where  N  is  an  adjustable  coefficient  (0  <  N  <  0.01 6)  and  x  is  distance  offshore.  Implicit  in  this  formula  is  that 
mixing  length,  and,  thus,  mixing  increases  with  distance  offshore. 

Kraus  and  Larson  (1991)  and  Smith  et  al.  (1993)  also  use  a  formula  for  eddy  viscosity  that  contains  an 
empirical  coefficient  A.  The  formula  assumes  that  mixing  is  primarily  a  function  of  orbital  wave  motions, 
i.e., 
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V,=A  umH\  (44) 

um  is  magnitude  of  the  horizontal  water  particle  velocity. 

Earle  (1991),  who  considered  three  formulas  for  eddy  viscosity,  accomplished  the  inclusion  of  lateral 
mixing  in  NSSM.  The  formula  given  by  Eq.  (25)  was  selected  because  the  dependence  on  horizontal  eddy 
viscosity  on  the  wave  breaking  dissipation  rate  is  intuitively  satisfying.  The  M  value  used  in  SURF  3.0  is  2. 

Changes  in  M  yield  changes  in  the  shape  of  the  V  profile  and  in  the  magnitude  of  the  highest  estimated 
value  of  V.  In  view  of  the  importance  of  M  in  producing  accurate  estimates  of  v,  this  value  is  re-examined  in 
this  report  using  a  greater  number  of  field  cases  than  used  in  previous  studies. 

5.0  LONGSHORE  CURRENT  MODEL  IMPROVEMENTS 

A  characteristic  problem  in  SURF  3.0  is  that  currents  are  generally  too  weak  in  the  main  portion  of  the 
surf  zone  but  often  unreasonably  high  in  very  shallow  water.  This  first  deficiency  is  demonstrated  using 
simple  statistical  analysis  of  model  errors.  By  showing  that  model  estimates  of  wave  height  are  largely 
correct,  it  is  inferred  that  energy  dissipation  processes  are  treated  adequately.  Therefore,  the  wave  momen¬ 
tum  input  is  correct,  and  the  deficiency  is  from  the  longshore  current  model  itself. 

Careful  troubleshooting  of  the  actual  model  code,  tests  of  the  algorithms  on  which  the  software  is  based, 
e.g.,  numerical  solutions  of  differential  equations,  and  review  of  the  theory,  revealed  no  apparent  underlying 
causes  for  these  model  deficiencies.  Simply  decreasing  cf  to  another  constant  value  does  not  produce  accu¬ 
rate  results  because  doing  so  does  not  effectively  change  the  shape  of  the  longshore  current  profile.  That  is, 
lowering  cf  will  increase  currents  over  most  of  the  surf  zone  but  at  the  expense  of  exacerbating  the  unreason¬ 
ably  high  currents  in  very  shallow  water.  Recent  measurements  of  c{  imply  that  a  variable  bottom  friction 
function  instead  of  fixed  constant  may  be  necessary.  But  no  existing  formulation  can  produce  consistently 
good  results. 

A  first  step  toward  an  improved  algorithm  is  to  recognize  that  the  maximum  longshore  current  is  often 
near  the  middle  of  the  surf  zone.  This  has  often  been  an  approximation  in  earlier  models  theoretically  appli¬ 
cable  to  planar  beaches  (e.g.,  Longuet-Higgins,  1970a,  1970b).  Lundberg  et  al.  (1997)  used  the  fact  that 
often  over  a  barred  beach  the  maximum  longshore  current  exists  near  mid  surf  zone.  This  phenomenon  is 
especially  evident  in  the  records  from  DELILAH,  in  which  the  current  maximum  is  often  over  the  trough, 
where  gradient  of  radiation  stresses  are  at  a  minimum. 

Figure  7  summarizes  the  maximum  longshore  current  locations  normalized  by  computed  surf  zone 
width  for  all  DELILAH  and  DUCK94  cases.  Cases  in  which  the  edge  of  the  surf  zone  is  estimated  to  extend 
seaward  of  the  model  starting  depth  are  not  used.  The  mean  ratios  from  DELILAH  and  DUCK94  of  the 
location  where  there  is  a  maximum  current  in  the  surf  zone  x(maxlVl)  to  the  width  of  the  surf  zone  Xh  are 
approximately  0.4  from  both  data  sets.  The  figure  shows  the  mean  value  of  the  ratio  (solid  line)  and  the  one 
standard  deviation  from  the  mean  (dotted  lines).  There  is  nearly  the  same  scatter  from  both  data  sets. 

The  use  of  the  model  to  estimate  surf  zone  width  is  based  on  the  recent  work  of  Mettlach  et  al.  (1999), 
who  showed  that  SURF96  accurately  estimates  surf  zone  width.  They  analyzed  658  video  images  to  obtain 
the  location  of  incipient  wave  breaking  and  then  compared  those  locations  to  model  output  to  derive  the 
optimal,  percent-breaking  threshold  from  the  Navy  model.  It  is  important  to  note  that  this  video  work  was 
done  using  SURF96  with  the  same  options  and  free  parameters  for  wave  transformation  and  energy  dissipa¬ 
tion  as  in  SURF  3.0:  rollers;  W(H)  in  the  form  of  Eq.(17);  y  =  0.42;  a  =  5°.  Surf  zone  width  is  estimated  to 
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DAY  (October  1994) 

Fig.  7  —  Time  series  of  relative  location  of  maximum  longshore  currents  within  the  estimated  surf 
zone.  Solid  horizontal  line  denotes  mean;  dotted  lines  denote  mean  plus  and  minus  one  standard 
deviation  from  the  mean. 


start  at  offshore  locations  where  10  percent  of  waves  are  estimated  to  be  breaking  or  broken.  From  the  video 
determinations  of  surf  zone  width,  using  a  10  percent  threshold,  the  model  produced  width  estimates  with  an 
expected  error  of  -2  m  and  one  standard  deviation  of  30  m  over  a  range  5  to  408  m. 

With  the  observation  that  a  current  maximum  exists  near  mid  surf  zone,  the  next  step  is  to  select  a 
functional  form  for  cythat  will  tend  to  produce  a  maximum  in  this  region.  Keeping  cf  a  constant,  relatively 
low  value  seaward  of  the  mid  surf  zone,  should  increase  the  maximum  current  estimate.  The  next  logical 
parameterization  of  cfis  to  steadily  increase  it  from  mid  surf  zone  to  the  edge  of  the  water.  This  is  accom¬ 
plished  by  making  cf  inversely  related  to  depth  h,  which  has  been  done  by  Smith  et  al.  (1993)  with  the 
Manning  friction  coefficient,  Eq.  (41),  and  by  Garcez  Faria  et  al.  (1998).  Increasing  friction  toward  the  water 
edge  effectively  reduces  currents  in  very  shallow  water. 

5.1  Development  of  a  Variable  Bottom  Friction  Coefficient  Function 

Based  on  the  facts  and  arguments  presented  in  the  previous  section,  the  following  depth-dependent 
function  for  c  was  developed: 


(  c* 

cf{x)  =  \  fhjXJ  2f 

f°{  Kx)  J 


\x>XJ  2 
\x<XJ2‘ 


(45) 
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Cp  is  a  constant  to  be  tuned  by  data,  and  Xb  is  surf  zone  width  computed  in  the  model.  The  function  cpc)  is 
based  on  Xb,  which  is  based  on  an  optimal  percent  wave  breaking  threshold  of  10  percent.  Figure  7  shows 
that  the  empirical  mean  value  of  is  approximately  0.40;  however,  a  value  of  0.50  is  used  in  Eq. 

(46). 

The  variation  of  ^inversely  with  depth  is  not  completely  arbitrary.  First,  there  is  preferential  sorting  of 
sediment  types  by  grain  size  in  the  nearshore,  with  larger  grain  sizes  in  the  intertidal  zone  and  finer  material 
offshore.  Howd  and  Birkemeier  (1987)  and  Birkemeier  (1991)  have  shown  that  at  Duck  beach  median  grain 
size  varies  from  0.70  mm  on  the  steep  foreshore,  to  0.20  mm  over  the  bar,  and  0.12  mm  offshore.  Second, 
Eq.  (45)  has  a  form  similar  to  the  effective  depth-dependent  functions  used  by  Smith  et  al.  (1993)  and  used 
by  Garcez  Faria  et  al.  (1998). 

The  typical  variation  of  the  proposed  clover  a  barred  beach  is  shown  in  Fig.  8,  which  gives  the  depth  and 
the  associated  function  c  for  the  2200  EST  10  October  1990  DELILAH  case.  The  edge  of  the  surf  zone  at 
Xb  =  128  m.  The  bottom  friction  cf{ x)  =  C/D=  0.003  is  constant  to  a  distance  of  64  m  with  depth  h-  1.6  m. 
Shoreward  of  mid  surf  zone  XJ2 ,  the  depth  increases  and  the  bottom  friction  decreases  to  a  minimum  of  cf= 
0.0024  at  x  =  41.2  m  and  h  =  2.1  m.  Bottom  friction  then  increases  linearly  near  the  edge  of  the  water  to  a 
maximum  of  cf  =  0.0141  at  x  =  3.4  m,  h  -  0.35  m.  The  model  automatically  terminates  at  the  most  shoreward 
grid  point  where  the  depth  is  greater  than  0.5  ft  (0.15  m).  Thus,  for  depths  greater  than  0.15  m  the  bottom 
friction  will  never  reach  more  than  about  ten  times  if  h(XJ2 ),  the  depth  at  mid  surf  zone,  is  1.5  m. 


Fig.  8  —  Profiles  of  depth  (top)  corrected  for  tide  stage  for  DELILAH  case  2200  EST,  10  October  1990,  and 
associated  function  ^(bottom)  using  C ^  =  0.003.  Surf  zone  width  for  this  case  is  estimated  by  model  to  be 
128  m. 
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5.2  Optimal  Cp  and  M 

The  basic  formulation  [Eq.  (46)]  yields  significant  improvements  in  the  current  if  the  proper  value  of  C)f1 
and  M  is  used.  Because  of  the  large  data  set  available,  various  combinations  of  these  two  parameters  are  used 
in  the  model  to  produce  estimates  for  each  of  the  DELILAH  and  DUCK94  cases,  which  are  then  compared 
to  measurements  from  which  model  skill  is  determined. 

Before  selecting  optimal  Cp  and  M,  it  is  useful  to  examine  the  effects  of  varying  Cp  and  M,  as  illustrated 
for  a  single  DELILAH  case  in  Fig.  9.  It  shows  several  current  profiles  from  surf  model  runs  for  which  or 
M  is  held  constant,  respectively.  The  increase  of  horizontal  eddy  viscosity  causes  more  mixing  and  produces 
a  less  peaked  longshore  current  profile.  The  mixing  depends  on  the  dissipation  rate,  therefore  it  has  the 
strongest  effect  near  the  bar  crest  where  dissipation  rate  is  often  highest.  The  parameter  Cfi]  has  a  strong 
effect  on  the  magnitude  of  the  current  throughout  the  whole  surf  zone. 


(a) 


(b) 


Fig.  9  —  Profiles  of  longshore  current  V  from  SURF  3.0  with  the  addition  of  Eq.  (45)— DELILAH,  10 
October  1 990, 2200  UTC.  Surf  zone  measurements  given  by  dots,  (a)  M  =  1 .0  for  all  profiles,  C  as  noted;  (b) 
Cfi  =  0.003  for  all  profiles,  M  as  noted 


Both  DELILAH  and  DUCK94  measurements  are  used  to  derive  the  optimal  values  of  and  M.  The  79 
DELILAH  cases  and  144  DUCK94  cases  were  selected  on  the  basis  of  data  availability.  Only  those  times  for 
which  wind,  wave,  depth,  tide,  and  validating  surf  zone  current  measurements  existed  were  used,  and  all 
such  cases  were  used.  Sixteen  sets  of  model  runs  were  made  using  different  combinations  of  Cfl  and  M 
parameters.  Model  skill  for  each  set  of  runs  is  summarized  in  Table  2  in  terms  of  various  skill  statistics.  The 
results  of  SURF  3.0  using  a  fixed  Cfll  are  listed  at  the  bottom  of  the  table.  This  case,  along  with  the  case  of 
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Table  2  —  Performance  of  longshore  current  estimates  when  using  denoted  C ^  in  the  variable  function  of 
cfx)  an  M.  Two  minimum  speeds  (0.001  m/s  and  0.1  m/s),  below  which  dtaa  are  not  included  in  computing 
the  statistics,  are  used  to  reduce  the  effect  of  large  number  of  data  points  at  low  current  speeds.  Bottom  row 
shows  the  results  from  the  old  model  with  constant  c r 
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C„  =  0.003  and  M  =  1 .0  are  in  bold  font  for  easy  comparison.  It  should  be  noted  that  large  numbers  of  low 
and  near-zero  current  speeds  can  contribute  high  relative  and  percentage  errors.  In  Table  2,  two  minimum 
speeds  (0.001  m/s  and  0. 1  m/s),  below  which  the  data  are  not  included,  are  used.  Under  the  higher  minimum 
speed,  the  total  number  of  points  drops  from  1922  to  1582.  Except  the  relative  errors,  most  of  the  statistics 
do  not  vary  much.  It  is  evident  that  decreasing  in  general  produces  better  statistics.  Among  several  best 
available  combinations,  =  0.003  and  M  =  1.0  were  selected  because  of  its  lower  £100and  standard  devia¬ 
tion  values.  Comparing  to  SURF  3.0,  the  e100  error  of  the  improved  model  has  reduced  to  17.5  from  52.3%. 
The  slope  of  the  linear  regression  curve  has  also  improved  from  0.48  to  0.89. 

From  the  preceding  analysis  of  model  errors,  summarized  in  Table  2,  the  optimal  values  C/0=  0.003  and 
M  =  1.0  are  used  for  the  cases  studies  presented  in  the  remainder  of  this  report.  The  value  of  =  0.003  is 
found  to  be  close  to  those  values  used  in  the  recent  literature. 

5.3  Scatter  Plots 

Figure  10  is  a  plot  of  1930  measurements  vs  improved  model  estimates  using  the  optimal  bottom  fric¬ 
tion  and  horizontal  mixing  coefficients.  (Skill  statistics  in  the  figure  are  defined  the  same  as  in  Table  2.) 
Compared  to  Fig.  4,  the  improved  model  achieves  significant  improvement  for  a  current  larger  than  0.5  m/s. 
The  improvement  is  also  evident  from  the  slope  of  the  linear  regression  line  changing  from  0.48  to  0.89. 


Cf=0 . 0030 ;  M=1.00 


Fig.  10  —  Scatter  diagram  of  measured  longshore  currents  vs  the  corresponding 
improved  model  estimates  of  current.  Dots  denote  DELILAH  measurements; 
circles  denote  DUCK94  measurements.  N  =  1930  observations.  Skill  statistics 
based  on  model  minus  measured  differences  of  set  from  combined  data  sets 
are:  linear  correlation  coefficient  R  =  0.89;  minimum  error  e  =  -1.06  m/s; 
mean  error  £  =  -0. 14  m/s;  maximum  error  e  =0.91  m/s;  standard  deviation 

mean  max 

of  the  errors  a  =  0.22  m/s.  The  solid  line  is  a  linear  regression  curve  with 
slope  m  =  0.89  and  y-intercept  =  -0. 14  m/s. 
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As  stated  before,  the  accuracy  of  the  model  in  estimating  the  maximum  current  in  the  surf  zone  is  very 
important  with  regard  to  modified  surf  index  (MSI)  computations,  for  which  only  the  maximum  current  in 
the  surf  zone  is  used.  The  maximum  longshore  current  from  the  two  data  sets  is  given  in  Fig.l  1,  which  can 
be  readily  compared  to  Fig.  5.  There  is  a  significant  improvement  in  estimates  of  V^;  the  root-mean-square 
error  decreases  from  0.32  m/s  to  0.27  m/s,  and  the  mean  fraction  error  decreases  from  -40  to  -30  percent.  As 
previously  mentioned,  the  statistical  improvement  should  be  a  lot  better  if  SURF  3.0  excludes  those  cases 
where  V  occurs  near  water  edge  rather  than  at  mid  surf  zone. 


Pig  ii  —  Scatter  diagram  of  maximum  longshore  current  measurement  within  surf  zone 
to  corresponding  improved  model-estimated  longshore  current.  Dots  denote  DELILAH 
measurements;  circles  denote  DUCK94  measurements.  N  =  223  observations.  Skill  statistics 
based  on  model  minus  measured  differences  of  set  from  combined  data  sets  are:  linear 
correlation  coefficient  R  =  0.899;  minimum  error  £mjn  =  -0.99  m/s;  mean  error  Zmean  =  -0.17 
m/s;  maximum  error  Ermx  =  0.41  m/s;  standard  deviation  of  the  errors  a=  0.27  m/s.  The 
solid  line  is  a  linear  regression  curve  with  slope  m  =  0.89  and  y-intercept  b  =  -0.15  m/s. 


6.0  ADDITIONAL  RESULTS 

In  this  section,  the  effectiveness  of  our  new  function  and  empirically  selected  and  M  is  shown,  using 
measurements  acquired  from  NSTS,  DELILAH,  and  DUCK94.  The  four  NSTS  cases  are  used  to  validate 
the  improved  model  because  no  data  were  used  in  selecting  Cp  and  M.  This  section  also  presents  several 
cases  studies  from  the  respective  DELILAH  and  DUCK94  data  sets. 
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6.1  Validation  Using  NSTS  Santa  Barbara  Data 

Results  from  four  cases  from  the  NSTS  at  Leadbetter  Beach,  Santa  Barbara,  California,  in  February 
1980  are  presented  next.  In  this,  and  the  two  following  sections,  model  results  from  SURF  3.0  with  and 
without  the  improved  bottom  friction  function  and  viscosity  parameter  are  examined.  The  improved  version 
was  run  using  =  0.003  and  M  =  1.0. 

Figures  12-15  give  cross-shore  profiles  of  depth,  root-mean-square  wave  height,  the  percentage  of  break¬ 
ing  waves,  and  longshore  current.  In  each  figure,  circles  denote  measurement  locations  or  measured  data.  In 
these  four  figures  and  in  similar  figures  in  Sections  6.3  and  6.4,  which  give  individual  cases  from  DELILAH 
and  DUCK94,  the  edge  of  the  surf  zone  is  denoted  as  a  dotted  vertical  line  in  the  third  panel.  This  is  the 
location  where  the  percentage  of  breaking  waves  is  estimated  to  reach  10  percent.  Relevant  offshore  wave 
parameters  used  for  model  initialization  are  noted  on  each  figure  beneath  the  top  panel. 

Thornton  and  Guza  (1986)  used  the  four  cases  in  this  section  to  validate  a  precursor  to  the  model  used 
in  SURF  3.0.  Earle  (1989)  used  three  of  the  four  cases  to  validate  the  first  version  of  the  NSSM  based  on 
Thornton  and  Guza  (1986).  These  are  ideal  cases.  First,  the  nearshore  depth  contours  during  the  experiment 
were  nearly  straight  and  parallel  to  the  beach.  Secondly,  the  offshore  waves  consisted  mainly  of  swells  with 
narrow  banded  frequency  and  direction. 


NSTS  Santa  Barbara  Feb.  03,  1980 


x  (m) 


Fig.  12  —  NSTS  planar  beach,  3  February  1980.  Profiles  of  surf  zone  parameters:  depth; 
root-mean- square  wave  height;  percent  breaking  or  broken  waves;  longshore  current,  standard 
model  (dotted)  and  improved  model  (solid).  Improved  model  estimates  based  on  C  = 
0.003  and  M  =  1.0. 
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NSTS  Santa  Barbara  Feb.  04,  1980 
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Fig.  13  —  NSTS  planar  beach,  4  February  1980.  Profiles  of  surf  zone  parameters:  depth; 
root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore  current,  standard 
model  (dotted)  and  improved  model  (solid).  Improved  model  estimates  based  on  - 
0.003  and  M  =  1.0. 


NSTS  Santa  Barbara  Feb.  05,  1980 


Fig.  14  —  NSTS  planar  beach,  5  February  1980.  Profiles  of  surf  zone  parameters:  depth; 
root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore  current, 
standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates  based  on 
C  =  0.003  and  M=  1.0. 
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Fig.  15  —  NSTS  planar  beach,  6  February  1980.  Profiles  of  surf  zone  parameters:  depth; 
root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore  current,  standard 
model  (dotted)  and  improved  model  (solid).  Improved  model  are  based  on  C_  =  0.003  and 
M=  1.0. 

Of  importance,  Santa  Barbara  data  have  not  been  used  in  model  calibration  in  choosing  the  bottom 
friction  and  eddy  viscosity,  so  it  is  indeed  encouraging  that  the  improved  model  produces  excellent  results 
compared  to  SURF  3.0.  In  all  four  cases,  the  improved  model  reaches  a  peak  that  is  closer  in  space  to  the 
observed  peak  than  is  the  peak  in  SURF  3.0.  The  shape  of  the  profile  from  the  improved  model  compares 
relatively  well  to  the  shape  of  the  measured  profile. 

6.2  DELILAH  Cases 

Figures  16-20  are  plots  for  five  DELILAH  cases  from  8  to  12  October  1990.  The  location  of  pressure 
gages  (+’s)  and  current  meters  (o’s)  are  denoted  in  the  first  panel  of  each  figure.  The  same  symbols  are  used 
to  denote  respective  measured  data  in  panels  2  and  4.  Tide,  wind  speed,  and  direction  are  also  listed  under 
the  first  panel.  Wave  parameters  from  the  directional  wave  spectra  acquired  at  an  8-m-deep  array  of  pressure 
gages  0.5  nmi  offshore  and  used  as  input  to  surf  model  are  lifted  under  the  second  panel.  Validating  wave 
height  and  longshore  current  measurements  are  averaged  over  the  same  time  duration  as  the  136-min  time 
duration  during  which  the  wave  spectra  were  acquired. 

Figures  16  and  17  represent  low  longshore  current  cases  with  no  obvious  bar  feature  on  the  beach 
profile.  The  new  model  shows  a  significant  improvement  in  longshore  current  prediction.  It  is  typical  that 
the  old  model,  as  in  Fig.  17,  shows  a  monotonic  increasing  longshore  current  toward  the  beach  with  a  high 
tide.  The  new  variable  bottom  friction  function  has  eliminated  this  incorrect  trend.  Figures  18  and  19  repre¬ 
sent  high  longshore  currents  over  a  barred  beach.  These  days  are  also  selected  because  the  wave  heights  are 
near  the  critical  value  for  some  amphibious  operations.  The  new  model  is  again  proven  significantly  better  in 
estimating  the  maximum  current.  It  also  eliminates  the  second  peak  near  the  beach.  Figure  20  shows  a  case 
of  low  current  over  a  barred  beach.  Compared  to  the  old  model,  the  location  of  the  maximum  current  is 
significantly  improved. 


(m/s)  PERCENT  H{m) 
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DELILAH — 9010081300 


Fig.  16  —  DELILAH  beach,  1300  EST,  8  October  1990.  Profiles  of  surf  zone  parameters: 
depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore  current, 
standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates  based  on 
C  =0.003  and  M  =1.0. 


DELILAH — 9010091000 
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Fig.  17  —  DELILAH  beach,  1000  EST,  9  October  1990.  Profiles  of  surf  zone  parameters: 
depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore  current, 
standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates  based  on 
C„  =  0.003  and  M=  1.0. 
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DELILAH— 90101016  00 


Fig.  18  —  DELILAH  barred  beach,  1600  EST,  10  October  1990.  Profiles  of  surf  zone 
parameters:  depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore 
current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates  based 
on  CL  -  0.003  and  M  =  1.0. 
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DELILAH — 9010111600 


Fig.  19  —  DELILAH  barred  beach,  1600  EST,  11  October  1990.  Profiles  of  surf  zone 
parameters:  depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves; 
longshore  current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model 
estimates  based  on  =  0.003  and  M  =  1.0. 
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Fig.  20  —  DELILAH  barred  beach,  1300  EST,  12  October  1990.  Profiles  of  surf  zone 
parameters:  depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore 
current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates 
based  on  CL  =  0.003  and  M  =  1.0. 
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For  all  cases,  including  Eq.  (46)  in  the  model  improves  the  overall  accuracy  of  the  estimated  current 
profile,  especially  the  maximum  current  in  the  surf  zone.  In  all  cases  however,  the  location  of  the  measured 
current  maximum  is  shoreward  of  that  estimated  by  the  model.  It  has  been  a  characteristic  of  the  DELILAH 
data  set  that  the  maximum  longshore  current  occurs  over  the  trough  shoreward  of  the  bar  crest.  But  the 
computed  maximum  occurs  over  the  crest  of  the  bar,  where  there  is  a  sharp  increase  in  wave  breaking  and 
energy  dissipation.  Much  research  by  the  writers  and  others  has  been  devoted  to  resolving  this  discrepancy. 
This  issue  is  additionally  discussed  later. 

6.3  DUCK94  Cases 

Figures  21-24  present  model  and  data  comparisons  selected  from  four  different  days  (October  10-13) 
of  the  DUCK94  experiment.  Waves  were  coming  from  northeast  for  October  10-12,  and  turned  to  be  from 
the  southeast  on  October  13.  Data  from  October  10-12  were  chosen  because  of  the  strong  longshore  current 
and  also  because  field  observations  of  bed  shear  stresses  are  available  (Garcez  Faria  et  al.,  1998).  October  1 3 
is  selected  to  provide  an  example  for  a  current  flowing  in  an  opposite  direction.  Perhaps  the  most  important 
difference  between  first  three  examples  and  the  DELILAH  cases  is  that  the  wave  forcing  is  from  the  north¬ 
east  rather  than  from  the  southeast.  Thus,  there  is  no  chance  that  current  forcing  is  perturbed  by  the  presence 
of  the  FRF  pier  and  the  deep  scour  hole  below  it,  which  is  located  south  of  the  surf  observation  area.  During 
all  four  days,  bar  features  exist  on  all  beach  profiles. 
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DUCK94  —  9410101600 


Fig.  21  —  DUCK94  barred  beach,  1600  EST,  10  October  1994.  Profiles  of  surf  zone 
parameters:  depth;  root-mean- square  wave  height;  percent  breaking  or  broken  waves; 
longshore  current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model 
estimates  based  on  Cm  =  0.003  and  M  =  1.0. 


DUCK94  —  94101116  00 


Fig.  22  —  DUCK94  barred  beach,  1600  EST,  11  October  1994.  Profiles  of  surf  zone 
parameters:  depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves; 
longshore  current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model 

estimates  based  on  C„  =  0.003  and  M-  1.0. 
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DUCK94—  9410111600 


Fig.  23  —  DUCK94  barred  beach,  1300  EST,  12  October  1994.  Profiles  of  surf  zone 
parameters:  depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore 
current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates 
based  on  C_  =  0.003  and  M-  1.0. 

p 


DUCK 94  —  9410131900 


Fig.  24  —  DUCK94  barred  beach,  1900  EST,  13  October  1994.  Profiles  of  surf  zone 
parameters:  depth;  root-mean-square  wave  height;  percent  breaking  or  broken  waves;  longshore 
current,  standard  model  (dotted)  and  improved  model  (solid).  Improved  model  estimates 
based  on  C  =  0.003  and  M-  1.0. 
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Figures  21  and  22  (October  10  and  11)  have  similar  wave  and  wind  conditions.  The  improved  model 
overpredicts  the  maximum  current  slightly,  but  the  unrealistic  second  peak  near  the  beach  from  the  old 
model  is  significantly  reduced.  Figure  23  (October  12)  represents  a  high  wave  event  with  significant  wave 
height  of  2.3  m  and  wind  speed  of  14  m/s.  Even  with  higher  waves,  the  longshore  current  is  smaller  than  on 
October  10-11  due  to  smaller  wave  angle  relative  to  beach  normal.  The  measured  current  peak  moves  further 
into  the  trough.  In  Fig.  24  (October  13),  the  wave  direction  is  coming  from  the  southeast  and  the  wind  has 
decreased  to  10  m/s.  Although  the  new  model  does  not  produce  a  maximum  longshore  current  as  high  as 
measured,  the  maximum  from  the  new  model  is  twice  that  of  the  old  model. 

In  general,  Figs.  21-24  show  that  the  improved  model  provides  better  longshore  current  estimates  than 
the  old  model.  Although  there  is  a  gradual  decrease  in  the  magnitude  of  the  computed  current  over  the  most 
shoreward  quarter  of  the  surf  zone,  measured  current  shows  much  sharper  decrease.  This  suggests  that,  for 
the  last  quarter  of  the  beach,  flow  resistance  begins  earlier  and  at  a  higher  magnitude  than  that  which  can  be 
controlled  by  the  improved  bottom  friction  function.  Other  flow  resistance  mechanisms  such  as  surf  zone 
turbulence,  i.e.,  the  turbulent  Reynolds  stress,  may  play  an  important  role.  The  turbulence  from  breaking 
waves  may  effectively  retard  mean  current  motion  due  to  vertical  shearing  stress  that  is  not  treated  in  one¬ 
dimensional  models.  It  should  be  pointed  out  that  the  ratio  of  turbulence  penetration  depth  to  local  depth  is 
relatively  large  at  very  shallow  depth,  making  this  omission  more  serious. 

7.0  DISCUSSION 

This  section  discusses  three  important  aspects  of  the  improved  longshore  current  model.  First,  the  dif¬ 
ference  between  horizontal  eddy  viscosity  values  used  in  our  one-dimensional  (ID)  and  another  quasi-3D 
model  is  discussed.  Next,  our  new  bottom  friction  function  is  compared  with  field  measurement.  Finally,  the 
location  of  the  longshore  current  maximum  is  discussed. 

7.1  Lateral  Eddy  Viscosity 

In  our  model,  the  turbulent  eddy  viscosity  is  based  on  Battles’  (1975)  formulation  as  shown  as  Eq.  (25). 
His  formulation  has  been  widely  used  by  nearshore  hydrodynamic  models.  The  empirical  constant  M  is 
generally  thought  to  be  of  order  one.  Based  on  laboratory  data,  Battjes  (1983)  showed  that  M  is  equal  to  0.8. 
A  wide  range  of  M  values  from  0.025  to  2  has  been  used  by  different  models.  The  optimal  value  we  derive  is 
1.0  whereas  SURF  3.0  uses  2.0.  In  a  recent  development  of  the  quasi-3D  SHORECIRC  model,  Sancho 
(1999)  used  an  M  value  of  0.1,  which  is  an  order  of  magnitude  smaller  than  ours.  There  are  at  least  two 
reasons  that  we  need  to  use  a  larger  value.  First,  we  need  to  compensate  for  the  limitation  of  our  ID  model 
in  not  including  the  nonlinear  interactions  between  cross-shore  and  longshore  currents.  This  mean  flow 
interaction  acts  as  an  additional  dispersive  mechanism  (Svendsen  and  Putrevu,  1994).  Second,  our  ID  model 
includes  only  lateral  or  horizontal  diffusion  through  v^but  not  the  vertical  diffusion,  i.e.,  vtz .  On  the  other 
hand,  SHORECIRC  considers  both  eddy  diffusions. 

7.2  Bottom  Friction  Function 

It  is  useful  to  compare  our  bottom  friction  function  with  field  measurements.  During  DUCK94,  current 
profiles  were  measured  using  a  specially  designed  sled.  Twenty -two  vertical  longshore  current  profiles  were 
taken  during  the  daylight  hours  of  October  10-12.  For  the  bottom  boundary  layer,  logarithmic  profiles  were 
fitted  to  the  data  based  on  a  least-squares  method.  cf  values  were  then  derived  from  the  slope  of  the  least- 
square  fit.  These  derived  cf  values  are  listed  in  Garcez  Faria  (1997). 
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In  Fig.  25,  daylight  hour  profiles  of  our  cf  and  the  derived  values  are  plotted.  The  cross-shore  distances 
given  in  Garcez  Faria  (1997)  are  converted  to  our  convention,  which  is  that  zero  is  the  edge  of  the  water 
based  on  the  depth  profile  and  tide  level.  We  do  not  know  the  exact  time  of  each  sled  run;  therefore,  mea¬ 
sured  revalues  are  plotted  as  an  error  bar  indicating  the  uncertainty  of  the  cross-shore  distance  due  to  the  tide 
stage. 


x(m) 

Fig,  25  —  DUCK94,  Comparison  of  cf  from  Eq.  (46)  and  measured  data  by  Garcez  Faria  (1997) 
(heavy  lines  give  the  cross-shore  range  of  the  measured  bottom  friction  coefficient;  vertical  lines 
denote  esimated  surf  zone  width  for  the  respective  case) 


It  is  encouraging  to  see  that  our  chosen  value  0.003  for  C '  over  the  more  seaward  half  of  the  surf  zone 
is  near  the  average  values  of  measured  bottom  friction  coefficeients.  This  value  is  not  inconsistent  to  the 
1986  Duck  measurements  of  Whitford  and  Thorton  (1996).  They  reported  that  mean  bottom  friction  coeffi¬ 
cient  are  0.004,  0.003,  and  0.001  for  offshore  the  bar,  on  top  and  immediately  in  front  of  the  bar,  and  in  the 
trough.  As  shown  in  Fig.  25,  our  friction  function  does  not  catch  the  rise  of  c,  soon  enough  shoreward  of  the 
current  peak  over  the  shoreward-half  of  the  surf  zone.  Also,  the  data  from  Garcez  Faria  show  a  fall  of  c, 
values  towards  shore,  but  no  current  profiles  farther  near  shore  were  taken  to  validate  our  assumption  that 
bottom  friction  increases  in  very  shallow  water.  The  lack  of  observations  at  such  depths  is  likely  due  to  the 
technical  difficulties  in  placing  sensors  in  such  very  shallow  water.  In  very  shallow  water,  the  penetration  of 
breaking  wave  occupies  most  of  the  water  column,  making  the  logarithmic  bottom  boundary  layer  hard  to 
define  and  measure.  Thus,  the  characteristics  of  c, in  very  shallow  water  are  not  well  known. 

There  are  reasons  that  justify  our  proposed  rise  in  c  near  the  water  edge.  First,  it  is  well  known  that  there 
is  a  general  increase  in  the  sand  size  from  seaward  to  the  water’s  edge,  especially  in  the  intertidal  zone.  The 
rise  in  our  cyin  very  shallow  water  accounts  for  this  trend  in  bed  roughness.  Second,  our  higher  values  of  cf 
toward  the  water’s  edge  partly  compensate  for  the  missing  vertical  mixing  of  our  ID  model  as  discussed  in 
Section  7.1. 
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There  are  other  possibilities  of  choosing  the  variable  in  Eq.  (46)  rather  than  the  depth,  e.g.,  the  percent¬ 
age  of  breaking  or  wave  height.  After  testing  several  parameters,  the  depth  turns  out  to  be  the  most  robust 
choice.  It  is  well  known  that  the  breaker  height  is  proportional  to  depth  in  the  saturated  surf  zone  due  to  the 
depth-induced  breaking.  Therefore,  the  ratio  between  depths  in  Eq.  (46)  is  equivalent  to  wave  height  ratio  in 
the  very  shallow  water.  This  may  be  the  major  reason  why  the  variable  friction  function  is  successful  in 
compensating  the  missing  vertical  mixing. 

7.3  Location  of  Longshore  Current  Maximum 

As  mentioned,  maximum  measured  longshore  current  often  occurs  in  the  trough  of  a  barred  beach.  On 
the  other  hand,  the  maximum  gradient  of  radiation  stress  from  our  model  and  others  is  always  located  over 
the  bar  crest.  The  shallow  depth  at  the  bar  crest  causes  an  increase  of  depth-induced  wave  breaking,  thus 
producing  a  local  current  maximum  near  same  location.  Many  research  efforts  have  investigated  this  dis¬ 
crepancy.  Using  a  large  wave  tank  and  a  barred  bottom,  Reniers  and  Battjes  (1997)  found  that  the  longshore 
current  peak  was  located  over  the  bar,  which  agrees  with  radiation  stress  theory.  They  concluded  that  the 
occurrence  of  the  maximum  longshore  current  velocities  in  the  trough,  as  observed  during  DELILAH,  is  not 
the  result  of  pure  local  wave  forcing.  However,  these  laboratory  results  must  be  framed  within  the  context  of 
length  scale  before  applying  them  to  field  scenarios.  Smith  et  al.  (1993)  suggested  that  the  spatial  lag  be¬ 
tween  production  and  dissipation  of  turbulent  kinetic  energy  causes  a  shift  in  the  momentum  decay.  In  other 
words,  the  wave  momentum  from  the  breaking  waves  is  temporarily  stored  as  turbulent  momentum  and  is 
not  immediately  used  to  generate  longshore  current.  They  solved  the  TKE  equation  and  incorporated  the 
results  into  a  longshore  current  model.  Using  DELILAH  data,  their  model  showed  only  small  improvements 
in  moving  the  computed  maximum  toward  the  trough.  In  recent  years,  quasi-3D  nearshore  wave-driven 
current  models  were  developed.  Garcez  Faria  et  al.  (1995)  investigated  the  effect  of  cross-shore  current  on 
the  location  of  the  maximum.  Using  SHORECIRC,  Sancho  et  al.  (1995)  studied  the  effects  of  longshore 
nonuniformities  on  longshore  currents.  None  of  these  models  had  complete  success  in  predicting  the  loca¬ 
tion  of  maximum  current.  More  research  needs  to  be  done. 

8.0  CONCLUSIONS 

In  this  report,  we  develop  a  variable  bottom  friction  coefficient  function  based  on  surf  zone  width  and 
depth.  It  is  important  to  note  that  our  function  is  not  dependent  on  information  about  bottom  material  and 
roughness,  which  are  often  difficult  to  obtain,  especially  in  denied  areas.  The  physical  reasons  for  the  valid¬ 
ity  of  the  function  have  been  covered  in  details  in  the  discussion  section.  In  short,  the  function  reflects  the 
increase  of  bottom  friction  due  to  the  sediment  sorting  and  compensates  for  the  lack  of  vertical  eddy  diffusivity 
in  ID  models.  The  latter  becomes  more  important  in  the  very  shallow  water  due  to  the  intensive  breaking. 

A  comprehensive  set  of  data  from  two  field  experiments  at  Duck,  North  Carolina,  is  used  to  develop  and 
calibrate  a  variable  bottom  friction  function  that  effectively  increases  the  skill  of  longshore  current  esti¬ 
mates.  The  new  bottom  friction  function,  along  with  a  newly  selected  eddy  viscosity  coefficient,  signifi¬ 
cantly  improves  the  magnitude  and  location  of  the  maximum  current,  and  greatly  reduces  the  unrealistic 
high  current  near  the  water  edge.  The  new  model  also  produces  lower  root-mean-square  errors  and  mean 
relative  errors.  Application  of  the  new  model  to  a  beach  at  Santa  Barbara,  California,  also  shows  significant 
improvements  over  the  old  model  with  the  fixed  bottom  friction  coefficient.  This  independent  validation 
proves  that  the  improved  model  and  our  approach  are  robust. 

The  improvements  described  by  no  means  imply  that  no  further  work  on  a  longshore  current  model  is 
needed.  Our  model  can  be  further  improved  in  its  prediction  of  the  location  of  the  maximum  current  on 
barred  beaches.  The  decay  of  longshore  current  in  the  most  shoreward  quarter  of  the  surf  zone  can  also  be 
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refined.  However,  there  is  a  limit  as  to  how  far  one  can  improve  ID  models,  which  are  understandably 
limited  in  the  representation  of  all  important  physical  processes  in  the  surf  zone.  In  a  few  years,  the  present 
ID  model  may  be  replaced  by  a  quasi-3D  model,  which  covers  nearshore  dynamics  more  completely.  An 
additional  advantage  of  a  3D  model  is  its  possibility  to  predict  rip  currents.  Until  then,  our  improved  model 
constitutes  significant  progress  in  longshore  current  prediction  and  will  be  used  to  upgrade  SURF  3.0.  It 
should  be  noted  that  we  have  demonstrated  the  usefulness  of  three  data  sets  acquired  during  intensive  field 
experiments.  These  data  sets  have  been  formatted  in  such  a  way  that  validation  of  future  surf  model  im¬ 
provements  can  be  done  efficiently.  Any  new  longshore  or  nearshore  current  model  needs  to  show  consis¬ 
tent,  accurate  performance  with  various  data  sets  before  it  can  be  adopted  for  operational  use. 
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Appendix  A 

VALIDATION  STATISTICS 

The  following  definitions  for  various  elementary  skill  statistics  are  used  throughout  this  report.  All  error 
statistics  are  based  on  the  difference  x  .  ,-x  ,  where  x  is  based  on  the  measurement.  For  A  observations 

model  true  true 

the  mean  of  the  simple  error  £  is  given  as 
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Physical  units  are  the  same  as  x.  and  x.  .  The  percentage  of  the  mean  fractional  error  or  100  times  the 
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mean  of  the  fractional  error  of  a  set  of  model  estimates  is  given  as 
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where  the  quantity  xis  a  threshold  value  chosen  to  prevent  values  of  infinity.  In  this  report,  values  of  x  are  a 
very  small  fraction  of  the  range  of  values.  The  mean  deviation  from  the  true  value  is  given  as 
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It  is  important  to  point  out  that  large  numbers  of  low  and  near-zero  values  can  contribute  to  high  mean 
fractional  errors.  It  is  often  necessary  to  exclude  some  low  values  to  avoid  distorting  the  results. 

The  maximum  and  minimum  errors  in  a  data  set  are  defined  as 

£min=min  (Xmodel-XtrJ9  (A4) 

^max  “  max(xmodel  ~  xtrue  )  •  (A5) 


The  functions  min( )  and  max( )  denote  the  minimum  and  maximum  errors  in  a  set,  respectively.  Note  that 
the  function  min( )  may  yield  a  value  that  has  an  absolute  value  greater  than  the  absolute  value  from  the 
corresponding  function  max( ). 


The  standard  deviation  o  is  the  sample  standard  deviation,  normalized  by  N  -1  as  given  by 
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Physical  units  are  the  same  as  x.mgdel  and  xjlrue. 
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The  linear  correlation  coefficient  R  is  defined  as 

N 

X  (%i  ~x)model{xi  ~  X\rue 

R  =  ^~ - ,  (A7) 

(N  1)0  model^t  rue 


where  the  sample  standard  deviations  in  the  denominator  are  defined  by  Eq.  (A8)  and  the  overbar  in  the 
numerator  (and  in  subsequent  definitions)  indicates  a  mean  value.  R  ranges  from  -1  to  1.  The  symbol  o., 
where  j  =  { model,  true},  is  not  explicitly  used  in  this  paper,  as  in  the  notation  Gmodel  and  Omodel  in  Eq.  (A6), 


N 


l  =  1 _ 


N- 1 


The  slope  m  of  the  least-squares  line  through  the  set  of  data  {x  ,  is  given  as 


(A8) 


m 


/V 

i^Xi,true  ^irue  )(^7,m  ode!  X model) 


-  X,ruef 


i=l 


(A9) 


The  ^-intercept  b  is  given  as 


b  —  x. 


m  odel  mxtrue  • 


(A10) 


Appendix  B 
NOTATION 


b  —  y  intercept  of  linear  regression  curve 

c  —  wave  phase  speed  (m/s) 

cd  —  air- water  drag  coefficient 

cf  —  water-seabed  friction  coefficient 

cg  —  wave  energy  speed  (m/s) 

c.  —  indexed  value  of  c  during  numerical  integration  (m/s) 

co  —  wave  phase  speed  used  for  model  initialization  (m/s) 
d  —  hydraulic  radius  of  open  channel  flow  (m) 

/  —  wave  frequency  (Hz) 

f  —  peak  wave  frequency  (Hz) 

g  —  acceleration  of  gravity  (m/s2) 

h  —  still  water  depth  (m) 

h.  —  indexed  value  of  h  during  numerical  integration  (m) 
k  —  wave  number  (radians/m) 

Ica  —  apparent  bed  roughness  length  scale  (m) 

k  —  indexed  value  of  k  during  numerical  integration  (1/m) 

m  —  slope  of  linear  regression  curve 

n  —  Manning  resistance  coefficient  (s/m3). 

p(H)  —  Rayleigh  wave  probability  density  function  (1/m) 

pb(H)  —  breaking  wave  probability  density  function  (1/m) 

r  —  roughness  parameter  (m) 

u  —  wave  orbital  speed  (m/s) 
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—  magnitude  of  the  horizontal  water  particle  velocity  (m/s) 

—  velocity  of  roller  in  kinetic  energy  equation  (m/s) 

—  cross-shore  distance  coordinate,  positive  seaward  (m) 

—  alongshore  distance  coordinate,  positive  left  of  x  (m) 

—  vertical  coordinate,  positive  up  (m) 

—  area  of  roller  (m2) 

—  simplified  notation  used  in  numerical  solution  of  longshore  current  representing  factor  of 
horizontal  mixing  (kg/s) 

—  simplified  notation  used  in  numerical  solution  of  longshore  current  representing  factor  of 
bottom  shear  stress  [kg/(m2  s)] 

—  constant  factor  in  the  function  cf 

—  simplified  notation  used  in  numerical  solution  of  longshore  current  (N/m2) 

—  simplified  notation  used  in  numerical  solution  of  longshore  current 

—  roller  energy  density  (J/m2) 

—  wave  energy  density  (J/m2) 

—  simplified  notation  used  in  numerical  solution  of  longshore  current  representing 

factor  contributing  to  longshore  current  (m/s) 

—  root-mean-square  wave  height  @(m) 

—  breaker  height  (m) 

—  root-mean-square  wave  height  (m) 

—  indexed  value  of  H  used  during  numerical  integration  (m) 

—  significant  wave  height  (m) 

—  wave  length  (m) 

—  coefficient  used  in  lateral  mixing  vf 

—  number  of  grid  points  in  numerical  solution  of  longshore  current 

—  linear  correlation  coefficient 
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wave  spectral  energy  density  function  [m2/(Hz-radians)] 
frequency-integrated  wave  spectral  energy  density  (m2/radians) 
cross-shore  component  of  cross-shore  radiation  stress  from  rollers  (J/m2) 
cross-shore  component  of  cross-shore  radiation  stress  from  waves  (J/m2) 
along-shore  component  of  cross-shore  radiation  stress  from  waves  (J/m2) 
wind  speed  (m/s) 
longshore  current  (m/s) 

indexed  value  of  longshore  current  used  in  numerical  integration  (m/s) 

maximum  longshore  current  in  the  surf  zone  (m/s) 

weighting  function  associated  with  Rayleigh  wave  distribution 

surf  zone  width  (m) 

wave  energy  dissipation  (W/m2) 

percent  fractional  error  (see  F  ) 

mean  of  the  absolute  value  of  the  fractional  error 

indexed  value  of  8  during  numerical  integration  (W/m2) 

maximum  error  in  model  estimate  compared  to  measurement 

mean  error  in  model  estimate  compared  to  measurement 

minimum  error  in  model  estimate  compared  to  measurement 

coefficient  of  0(1)  used  with  Thornton  and  Guza  wave  dissipation  model 

wave  setup  (m) 

lateral  mixing  eddy  viscosity  used  in  SURF  3.0  (m2/s) 

density  of  seawater  (kg/m3) 

density  of  air  (kg/m3) 

angle  of  the  wave-roller  stress  vector 

standard  deviation 
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—  minimum  threshold  value  in  denominator  when  calculating  fractional  error 

—  TKE  production  at  wave-roller  interface  (J/m2) 

—  longshore  component  of  mean  time-averaged  bottom  shear  stress  (N/m2) 

—  wave  stress  in  the  longshore  direction  (N/m2) 

—  wind  stress  in  the  longshore  direction  (N/m2) 

—  4  tan1(l)  =  3.14159265. . . 

—  wave  angle  (direction  from;  positive  counterclockwise  from  shore  normal;  radians) 

—  average  wave  direction 

—  indexed  value  of  wave  angle  (direction  from;  positive  counter-clockwise  from  shore 
normal;  radians) 

—  wave  angle  used  for  model  initialization  (direction  from;  positive  counter-clockwise  from 
shore  normal;  radians) 

—  peak  wave  angle 

—  wind  direction  angle  (direction  from;  positive  counter-clockwise  from  shore  normal; 
radians) 

—  angular  frequency  (Hz) 

—  maximum  water  particle  excursion  amplitude  of  fluid  motion  at  the  bottom  (m) 

—  tide  level  (m) 

—  empirical  coefficient  used  in  Kraus  and  Larson  formula  for  eddy  viscosity 


